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Abstract 

Consider an arbitrary complex- valued, twice continuously differentiable, nonvanishing 
function ip defined on a finite segment [a, b] C HL Let us introduce an infinite system of 
functions constructed in the following way. Each subsequent function is a primitive of the 
preceding one multiplied or divided by <p alternately. The obtained system of functions 
is a generalization of the system of powers {{x — xo) k } kQ . We study its completeness 
as well as the completeness of its subsets in different functional spaces. This system of 
recursive integrals results to be closely related to so-called L-bases arising in the theory 
of transmutation operators for linear ordinary differential equations. 

Besides the results on the completeness of the system of recursive integrals we show 
a deep analogy between the expansions in terms of the recursive integrals and Taylor 
expansions. We prove a generalization of the Taylor theorem with the Lagrange form of 
the remainder term and find an explicit formula for transforming a generalized Taylor ex- 
pansion of a function in terms of the recursive integrals into a usual Taylor expansion. As 
a direct corollary of the formula we obtain the following new result concerning solutions 
of the Sturm-Liouville equation. Given a regular nonvanishing complex valued solution 
yo of the equation y" + q(x)y = 0, x € (a, 6), assume that it is n times differentiable at a 
point xq € [a, b] . We present explicit formulas for calculating the first n derivatives at xo 
for any solution of the equation u" + q(x)u = Xu. That is, an explicit map transforming 
the Taylor expansion of yo into the Taylor expansion of u is constructed. 

AMS subject classification: 34B24; 41A30; 42A65 

Keywords: complete system of functions; Sturm-Liouville problem 



1 Introduction 

In the recent work [16] by means of pseudoanalytic function theory [2] , [17] a representation 
for solutions of Sturm-Liouville equations in terms of spectral parameter power series was 
obtained. For a proof not requiring pseudoanalytic functions we refer to [T7] and [19]. The 
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representation turned to be an appropriate tool for solving different Sturm-Liouville and 
related problems |31,0,0,ini,II51,[I31,II3,II51,IISl,|2Dlaiid[3I]. As is well known (see, 
e.g., [23]) under certain regularity conditions a solution to the initial value problem for the 
Sturm-Liouville equation is an analytic function of the spectral parameter and hence admits 
a Taylor series expansion in powers of the spectral parameter. In fact in [TB] the recursive 
formulas for calculating the Taylor coefficients in that expansion were proposed. The Taylor 
coefficients are naturally functions of the independent variable. Considered as an infinite 
system of functions they form a so-called L-basis related to a linear differential operator 
L. The theory of L-bases and L-analytic functions was developed in [9] (unfortunately that 
interesting and important book has not been translated into English) in relation with the 
concept of the operators of transmutation (see, e.g., pQ, [4], [26] and the recent review [25]). 
also called in Russian bibliography operators of transformation (see, e.g., [22J, [24 J). 

In [9] it was shown that to every (regular) linear ordinary differential operator one can 
associate a linear space spanned on an L-basis, and the equivalence of those linear spaces 
corresponding to operators of the same order was proved. Then the operator of transmutation 
can be regarded as an operation which transforms functions from one such linear space 
corresponding to a certain operator L to functions from another linear space corresponding 
to another operator M, and the transformation consists in substituting the L-basis with the 
M -basis keeping the same coefficients in the expansion. In the recent work [18] it was shown 
that the situation with the linear space generated by the L-basis is in a sense simpler and 
more natural. Namely, the L-basis corresponding to a regular Sturm-Liouville operator on a 
finite interval (a, b) is complete in L2(a, b). Moreover, in the present work we prove that in 
the space of piecewise continuously differentiable functions on [a, b] the L-basis is complete 
with respect to the maximum norm (see theorem [9]). We show a deep analogy of the L-basis 
with the system {(x — xo) n }^L Q of powers of the independent variable. The formulas for 
calculating coefficients in the expansion of a given function in terms of those generalized 
powers resulted to be direct generalizations of the formulas for Taylor coefficients. We prove 
a generalization of the Taylor theorem with the Lagrange form of the remainder term. 

A natural question then is about a relation between the generalized Taylor coefficients of 
a function and its ordinary Taylor coefficients. This question has not been studied previously. 
We find this relation and use it to obtain the main result of this paper, the representation for 
ordinary Taylor coefficients of any solution of the Sturm-Liouville equation u" + q{x)u = Xu 
in terms of the Taylor coefficients of a particular solution uq of the equation y" + q(x)y = 0. 
More precisely we obtain the Taylor coefficients of the quotient u/yo. It is interesting to 
notice that the form of the map transforming one set of Taylor coefficients into the other is 
independent of the particular form of yo. The paper ends with examples of application of 
the presented result. 

2 Definition and example 

Let / G C 2 (a, b) n C 1 [a, b] be a complex valued function and f(x) ^ for any x S [a, b]. The 
interval (a, b) is supposed to be finite. Let us consider the following auxiliary functions 




(1) 
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X 

X<ri(x) = nf X^ l \s) (/ 2 ( S )) ( " ir_1 ds, (2) 

x 

X 

xW(x) = nf X^is) {f(s)) { - ir ds, (3) 

so 

where x$ is an arbitrary fixed point in [a, b]. We introduce the infinite systems of functions 

{fn}n=l and {9n}n=l ( 4 ) 

defined by the relations 

f n = fX^ n ~^ and g n = fX( 2n -V. (5) 

In |18j it was shown that both systems of functions in (JH) are complete in L 2 {a, b) when the 
point xq coincides with one of the end points of the interval, and if xq is an interior point of 
the interval, then the completeness of the union of two systems {/ n }^Li U {g n }^Li can be 
guaranteed. 

The question about the completeness of the systems (|3J) is natural due to the following 
observation. 

Example 1 Let f = 1, a = 0, 6=1. Then it is easy to see that choosing xq = we have 

h(x) = l, f 2 (x) = x 2 , h{x) = x\... (6) 

and 

gi(x)=x, g 2 {x) = x 3 , . . . . (7) 

As is well known due to the Miintz theorem (see, e.g., J21 p. 270]) both systems of polynomials 
are complete in L 2 (0, 1). The systems Op represent a direct generalization of the systems of 
polynomials (0|) and ^ if instead of f = 1 an arbitrary sufficiently smooth and nonvanishing 
function is chosen. 

The completeness of the system ([6]) can be regarded as a corollary of the completeness 
of the system of eigenfunctions of the Sturm-Liouville problem 

u" + Xu = 0, 0<x<l, (8) 

u'(0) = u'{\) = 0. (9) 

Indeed, the eigenfunctions of this regular Sturm-Liouville problem have the form u n (x) = 
cosnirx, n = 0, 1, 2, . . ., and as each of them admits a uniformly convergent Taylor expansion 
in even powers of x, the system Q is also complete. 

Analogously, considering (J8|) with the boundary conditions u(0) = u(l) = we arrive at 
the complete system of eigenfunctions v n (x) = sinnirx, n = 1, 2, ... . Each of them admits 
a uniformly convergent Taylor expansion in odd powers of x and hence the system ([7J is 
complete in £2(0, 1) as well. 
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3 Completeness of the systems of recursive integrals 



We will need the following theorem from [16] (for additional details see [T7] and [19]) estab- 
lishing the relation of the systems of functions (|3J) to Sturm-Liouville equations. 

Theorem 2 ]16$ Let q be a continuous complex valued function of an independent real vari- 
able x G [a, b], X be an arbitrary complex number. Suppose there exists a solution f of the 
equation 

f" + qf = (10) 

on (a, 6) such that f £ C 2 (a,b) together with 1/f are bounded on [a,b]. Then the general 
solution of the equation u" + qu = Xu on (a, b) has the form 

U = C\U\ + C2U 2 

where c\ and c 2 are arbitrary complex constants, 

\k 00 \k 



— A/ — 



with X( n ) and X^> being defined by (GP-(E]) and both series converge uniformly on [a,b]. 

Remark 3 It is easy to see that by definition the solutions u\ and U2 satisfy the following 
initial conditions 

ui(x Q ) = f(x ), u[(xo) = f'(x ), (12) 
u 2 (x ) = 0, ^(xo) = l/f(x ). (13) 

Remark 4 It is worth mentioning that in the regular case the existence and construction of 
the required f presents no difficulty. Let q be real valued and continuous on [a, b] . Then (ffl\) 
possesses two linearly independent regular solutions v\ and v 2 whose zeros alternate. Thus 
one may choose f = v% + iv 2 . 

Theorem 5 fl8\j Let (a, b) be a finite interval and f 6 C 2 (a, b) Pi C 1 [a, b] be a complex valued 
function such that f(x) ^ for any x £ [a, b]. Then both systems of functions defined by 
the relations |5]) and (Op-® with xq = a are complete in L 2 (a,b). 

Proof. Under the conditions of the theorem, on the interval (a, b) the function / is a regular 
solution of the equation (fit)]) with q := —f"/f being a continuous complex valued function. 
Let us consider the equation 

u" + qu = Xu (14) 

with the boundary conditions 

u(a) = u(b) = 0. (15) 

It is known (see \24\ p. 36]) that the system of all eigenfunctions and generalized eigenfunc- 
tions of this problem is complete in L 2 (a,b). Due to theorem [2] and remark [3] if X n is an 
eigenvalue of the problem (fl4"j) . (fT5l) then the corresponding eigenfunction up to a constant 
factor must coincide with u 2 from (jlip where A = A n , i.e., it admits a uniformly convergent 



4 



series expansion in terms of the system of functions g n from ([5]). Moreover, if the multiplic- 
ity of A n is greater than 1 then the corresponding generalized eigenfunctions are obtained 
differentiating the eigenfunction with respect to the spectral parameter A (see |244 p. 27]). 
The representation (jlip of U2 shows us that the result of this operation will be again a series 
in terms of the functions g n . Thus, the eigenfunctions and the generalized eigenfunctions 
of the problem (j!4|) . f)15|) can be represented as uniformly convergent series in terms of the 
system of functions {fltal^Li- Then due to the Lauricella theorem about the transitivity of 
the property of completeness (see, e.g., JSJ p. 264]) we obtain that the system {<7n}^Li is 
complete in £2(0, b). 

In a similar way the completeness of {/n}^Li is proved by considering the Sturm-Liouville 
problem for equation (|14p with the boundary conditions 



All eigenfunctions of this problem coincide with u\ from (jlip where A = A n (this is due 
to remark [3]) and generalized eigenfunctions are obtained from the eigenfunctions by dif- 
ferentiation with respect to A. Thus, by analogy with the previously considered case all 
the eigenfunctions and generalized eigenfunctions of this Sturm-Liouville problem are repre- 
sented as uniformly convergent series in terms of f n , and by the Lauricella theorem {f n }^Li 
is complete in £2(0, b). ■ 

In this theorem we assumed that the point xq coincided with one of the end points of the 
interval (a, b). When xq is an interior point of the interval in general the systems {/ n }^Li 
and {gn}^ = i separately are not complete. It is easy to see on the considered above example 
[I] that if a = —1 and all other values remain unchanged then the system ([6]) is not complete 
anymore in L2(a,b) because ([6]) contains even functions only. Nevertheless considering the 
combined system of functions {f n }^=i^{9n}^Li we obtain a complete system. The following 
theorem establishes that this remains true in a much more general situation. 

Theorem 6 fl8\/ Let (a, b) be a finite interval and f £ C 2 (a, b) n C 1 [a, b] be a complex valued 
function such that f{x) 7^ for any x £ [a, b]. Then the system of functions {/n}^LiU{<7„}^ ! L 1 
defined by the relations (GJj and (CP-(0j with xq being an arbitrary point of the interval [a,b] 
are complete in £2(0, b). 

Proof. As in the proof of the previous theorem let us consider equation (|14p where q := 
—f" I f 1 for example, with the boundary conditions (|15p . Due to theorem[2]any eigenfunction 
as well as any generalized eigenfunction of this problem can be represented as a uniformly 
convergent series in terms of the functions f n and g n . From the completeness of the system 
of eigenfunctions and generalized eigenfunctions in L2(a,b) and by means of the Lauricella 
theorem we obtain the result. ■ 

Remark 7 If instead of the "seed" function f one considers the function 1/f then, as it is 
easy to see the auxiliary functions X^ and change their respective roles, and in the 

same way as was done above one can prove the completeness of the systems 



in £2(0, b) when xq coincides with a or b and the completeness of the union of these two 
systems in £2(0, b) when xq is an arbitrary point of the interval [a, b\. 



f'(a)u(a) - f(a)u'(a) 



u{b) = 0. 




and 
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Remark 8 Under the considered conditions on the function f it is easy to prove the com- 
pleteness of the systems {X^™" 1 ))}^, {X^ 1 )}^ p {^ (2n - 1} }" =1 , {X^ n ~^ }~ 1 when 
xq coincides with a orb and the completeness of the systems |X^ 2 ( n_1 ^ }^ =1 U {X^ 2n ~ 1 ^}'^' =1 
and {^ (2n_1) }~ =1 U {^ (2(n_1)) }^ = l when x is an arbitrary point of the interval [a, b] . This 
is based on the previous theorems and on the observation that if{h n }™ =1 is a complete system 
in L,2(a,b) then taking a continuous in [a,b] weight function p, such that l/|p| is separated 
from zero one has that {phn}^ =1 is complete in L2(a,b) as well. For the proof see J7| Sect. 

4.7]. 

Let us introduce the system of functions {^fcj^o defined as follows 

f f(x)X^){x), k odd, 
Vkix) = { f(x)X^(x), k even, (16) 

where the definition of X^ and X^ is given by ([I])-© with xq being an arbitrary point of 
the interval [a, b]. We are interested in the completeness of this system in the space of piece- 
wise differentiable functions with respect to the maximum norm and in the corresponding 
series expansions. 

The system of eigenfunctions of the Sturm-Liouville problem under certain regularity 
conditions is complete not only in the sense of the L2-norm but also in the uniform conver- 
gence topology. The approach used in the proof of the previous two theorems and based on 
the completeness results for the system of eigenfunctions and generalized eigenfunctions of 
the Sturm-Liouville problem as well as on the Lauricella theorem about the transitivity of 
the property of completeness gives us the following result. 

Theorem 9 Let f satisfy the conditions of theorem^ and {^fcj^o ^ e ^ e s D s ^ m °f func- 
tions defined by U6\) with xq being an arbitrary point of the interval [a, b] . Then for any 
complex valued piecewise continuously differentiable function h defined on [a, b] and for 
any e > there exists such N £ N and such complex numbers at, k = 0, 1, ... N that 

N 

' < e. 



max z€[a,&] 



fc=0 



Proof. Take a piecewise continuously differentiable function h. It satisfies certain boundary 
conditions of the form 

dh(a) + c 2 h'(a) = and c 3 h(b) + c A h'(b) = (17) 

where c\ + c| ^ and c| + c\ ^ 0. Let {wfel^o ' 3e a system of all eigenfunctions and 
generalized eigenfunctions of the Sturm-Liouville problem for equation (|14p with the non- 
degenerate boundary conditions (|17|) . This system is complete in the linear space of piecewise 
continuously differentiable functions with respect to the maximum norm (see [24|, Chapter 
1]). As was previously shown every admits a uniformly convergent series expansion in 
terms of the functions (p^. Thus again, the result is a corollary of the Lauricella theorem. ■ 

Remark 10 In what follows instead of the system of functions {<Pk}kLo it will be slightly 
more convenient to consider the system {V'fc = "y-} fc _ which from the previous theorem and 
due to the boundedness of\f\ and l/\f\ is also complete with respect to the maximum norm 
in the space of piecewise continuously differentiable functions defined on [a, b] . 



6 



4 Generalized Taylor expansions 



In the previous section we showed that the system of recursive integrals 

f Xto(x), k odd, 

y>*W-\ x(k) {xh keven 

is complete in Li2(a, b) as well as with respect to the maximum norm in the space of piecewise 
differentiable functions defined on [a, b\. The system of recursive integrals generalizes the 
system of powers of the variable x. In this section we obtain generalizations of several 
classical results concerning the Taylor expansions. We start with a generalization of Taylor's 
formula. 

Let us assume that the functions / and h on a certain segment [a, b] possess the derivatives 
of all orders up to the order n and that f(x) 7^ 0, when x £ [a, b]. Then in [a, b] the following 
generalized derivatives are defined 

7o0)0) = h(x), 

lk{h){x) = (f\x)) { - 1)k " {7k-i(h))'(x), (18) 
for k = 1, 2, . . . , n. Let us consider a function of the form 

n 

Pn(x) = ^2a k Mx)- (19) 

fc=0 

Similarly to the fact that the coefficients of a polynomial £^ =0 afc(x — xo) k can be expressed 
through its value and the values of its derivatives at the point xo we obtain that the coeffi- 
cients «fc in (fT9|) can be expressed through the value of P n and the values of its generalized 
derivatives at the point xq. Indeed, a simple calculation gives us the following result 

7k(Pn)(x ) , . 

k = — ifei — ■ ( > 

Functions of the form (119p will be called generalized polynomials of order n. Notice that 
application of the generalized derivatives ()18|) to a generalized polynomial does not require 
the smoothness of /,- formula (I20p is true for any continuous and nonvanishing /. 

Now let us consider a function h possessing at the point xq the derivatives of all orders 
up to the order n and suppose that the same is true for the function / which additionally 
is different from zero at Xq. In relation with the function h we introduce a generalized 
polynomial of the form f)19|) where the coefficients are defined by the equality 

7k(h)(x ) 

k = kl • ( ' 

According to the previous observation this generalized polynomial together with its general- 
ized derivatives (up to the order n) at the point xq possess the same values as the function 
h and its derivatives, 7fc(-Pn)(£o) = Jk(h)(xo), k = 0, 1, ... n. We are interested in estimat- 
ing the difference between P n (x) and h(x) for x ^ xq. In the next theorem we obtain a 
generalization of the theorem on the Taylor remainder term in the Lagrange form. 
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Theorem 11 ( Generalized Taylor theorem with the Lagrange form of the remainder term ) 
Let {f,h} C C n+1 [xo,b] and f(x) 7^ there. Then for any x 6 [a;o,&] there exists a number 
c between x$ and x such that 

,/ \ lk{h)(x Q ) 7 w+ i(/t)(c) 

fc=0 ' 

Proof. Consider the difference i? n = h — P n and the function V'n+i- For simplification of 
the notation let us skip the subindices, R = R n and ip = ip n +i- We have that 7fc(i?)(xo) = 
7fc(V0( x o) = 0, A; = 0, 1, ... n. We may therefore apply Cauchy's mean value theorem to the 
functions R and 

R(x) = R(x) - R(x ) = R f ( Xl ) = fjx^R'jx,) = 7 i(i?)(ai) 
ip(x) ip{x)-ip(x ) ip'(xi) f 2 {xi) tp'(xi) 7i(V>)(a;i) 

where xq < x\ < x < b. Another application of Cauchy's mean value theorem gives us the 
equalities 

7i(it!)(xi) _ 7i(/2)(xi)- 7 i(/2)(xo) _ 7i(i2)(s 2 ) _ r 2 (x 2 ) 71(^2) _ 72(^)^2) 



7i(V?)(xi) 7i(^)(a;i) - 7i(^)(x ) 7iW( a; 2) / 2 (^2) 7iW(^2) 72(V0( :c 2) 
where xq < x% < x\. Continuing this procedure we obtain 

j n {R){x n ) = -/ n (R)(x n ) - 7 w ( J R)(x ) = j n (R)(x n+ i) 

1n{i>)(x n ) n/n(lp)(x n ) ~ 7n(^)(«o) TnW^n+l) 

_ /(- 1 ) n+l2 (x n+1 ) 7 ;(fi)(x n+1 ) _ 7 n+l(^)(x ra+ l) 
/(~D" +l2 (x n+1 ) 7;W(^n+l) 7n+l(^)(x n+ l) 
where xo < x n+ \ < x n < x n _i < ■ ■ ■ < x < b. Consequently, 

R(x) _ jn+i{R)(x n +i) 
ip(x) -y n+1 (if>)(x n+1 ) ' 

Moreover, by the definition of the functions ipk and of the generalized derivatives we have 
that j n+ i(R)(x) = j n+1 (h)(x) and j n+1 (ip)(x) = (n + 1)!. Then 

where c = x n+ i. ■ 

Obviously, the classical Taylor theorem with the Lagrange form of the remainder term is 
a special case of theorem 1111 when / = 1. 

Definition 12 A functional series of the form 

k=0 

will be called the generalized Taylor series of the function h. 
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Theorem [2] gives us an important example of uniformly convergent generalized Taylor 
series. Indeed, we obtain that the quotients of the linearly independent solutions u\, u<i of 
equation (|14p and the particular solution / of ([8]) have the form 



Ul 

f 



— A^n and = V 

n=0 J n=0 



1 + (_!)»+! 

2n! 2 ^ 



5 A relation between the generalized and the classical Taylor 
expansions 

In this section we establish relations between the generalized Taylor coefficients of a suffi- 
ciently smooth function and its usual Taylor coefficients. 

Theorem 13 Let {f,h} C C n [a,b] and f(x) / 0, x G [a,b] and ip = f 2 . The following 
relation between the ordinary derivatives and generalized derivatives of the function h at the 
point xq G [a, b] are valid 



( h(x ) \ 




h'(x ) 




h"(x ) 




h"'(x ) 




\ hW(x ) ) 
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l/ip(x ) 

a 2 ,i(x ) 1 

a 3j i(x ) a 3i2 (x ) 1/^(^0; 








!+(-!)" 



\ / 7 o(M(xo) \ 
li(h)(x ) 
l2(h)(x ) 
l3(h)(x ) 



"(*o) J 



where the functions a n ^ m are calculated following the recursive procedure: 

— ' _i_ (-!) m 

"n,m — ^n— 1 m ' V Q"n~ l,m— 1- 



\ 7n0)0o) J 

(22) 
(23) 



Proof. The first two rows in (|22p we obtain directly from the definition of the generalized 
derivatives. To obtain h" we proceed in the following way 



From (1181) we have 



{ lk {h))' = <p { - X) lk+i(h). (24) 

In particular, (71 (/t))' = <pj2(h). Thus, h" = 72(h) + 71(h) (and hence a 2 ,i = \ )■ 
In general, assume that we have calculated the row corresponding to h^ where k is even, 

h [k] = 7k(h) + ajfc,fe-i7fe-l(^) H r- afc,i7i(^)- 

Differentiating this equality we obtain 

/i [fc+1] = (7fc(M)' + o*,fc-i (7fc-i(M)' + • • • + a k ,i (11(h))' + 4 fc _ l7fe _i(/i) + • • • + a' k l7 i(/i). 
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Replacing the derivatives of Jj(h), j = 1, . . . , k with y>( -1 ) J jj + i(h) we arrive at (|23|) . For 
an odd A; the reasoning is analogous. ■ 

The transformation matrix in (|22p will be denoted by A n . 

Meanwhile the recursive procedure (|23p allows one to calculate the first few rows of the 
transformation matrix A n the next statement gives us a general formula for its elements. 

Proposition 14 The element a n ^ m with 1 < n < N and 2 < m < n of the matrix An, 
N = 2,3,... in {HP has the form 



0-n,rn — ^ ( , j (— j Ofc,m-l 
k=m-l ^ ' ^ ' 

where (p represents the binomial coefficients, 

e (^)^ s (v) (?r^ 

fc 1= m-l V 7 fc 2 =m-2 v * / \r/ 

fc m _ 3 -l . . . . [fc m _ 3 -l-fc m _ 2 ] fc m-2-l /, ^ /-i\[fem-l-l] 



«m-2 



fcm-2=2 

iu/ien m is even and 

<~- e e (v)(?r^- 

fc 1=m -l v 1 7 k 2 =m-2 V 7 XYV 

^-a- 1 /, ,\ fcm-a-1 /, 1N / 1 v [fc m _ 2 -l-fe m _i] 



e rr: K"- 1 -^ e t; u 

rirn-2 J t _i \ "W— 1 J X^P J 

fcm—l — J- 



^[fcm-l-l] 



u;/ien m is ode? fin 6oi/i cases, 2 < m < k). 
Proof. Let us notice that 

ti = -71 (25) 

<P 

(we will write 7^ understanding jk(h)) and recall the formula for the n-th derivative of a 
product of two functions 



Applying this formula in order to evaluate the (n — l)-th derivative of the equality (|25p we 
obtain 

n-l , 1 \ / 1 \ [n— 1— k] 

* )(^ ^ n> ~ h (26 » 
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Considering j[,k>lwe obtain 




Continuing this procedure we have 

fifcl 

7i = h,i72 + b k ,273 + frfc,374 H 1- h,k^k+i 



where 
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and in general bk t m are defined by the equalities from the statement of the proposition. 
From (|26p we obtain for n > 1 



where 



7i 



and hence 



-. \ [n-ll n-1 / i\ /i \ [n-l-fcl 
a n ,i7i + an,272 + a n ,373 + a™,474 H h a njn _i7„_i + a ni „7 n , 



! \ [«-l] 
ra— 1 / — 1\ / 1 \ [ n ~ 1—S; ] 

n— 1 / — l\ / 1 \ I n_1_fc ] 



n— 1 / — lA / 1 \ [ n_1_fe ] 

«re,n-l = / , ( , I — ^fc,n-2 



n— 1 / i \ / i \ [n— 1— fel 



k=n 



l n,m — 

k=m— 1 
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6 Taylor coefficients of solutions of the Sturm-Liouville equa- 
tion 



In the present section we obtain an interesting corollary of theorem [13j Given a nonvanishing 
(complex valued) solution of equation (fTOj) possessing n derivatives at a certain point, for 
any solution of (|14p (for any value of the spectral parameter A) we are able to calculate the 
exact values of its n derivatives at the same point without any integration. Moreover, each 
derivative is a polynomial of a certain order with respect to A. 

Theorem 15 Let (a,b) be a finite interval and f 6 C 2 (a,b) nC 1 [a, b] be a complex valued 
solution of n~U\) such that f(x) ^ for any x 6 [a, b] and at the point xq £ [a,b] there exist 
the derivatives of f up to the n-th order. Then the linearly independent solutions u\ and 
u 2 of equation (TJ\ ) satisfying the initial conditions U^l and i TO)) respectively possess at the 
point xq the derivatives up to the n-th order which can be calculated according to the following 
relations 



and 



( ui(x )/f(x ) \ 

(Ui/f Y (x ) 

(ui/f)" (xo) 
Oi//)"' ( x o) 
( Ul /f) IV (x ) 

\ («i//) W (*o) J 



( u 2 (x )/f(x ) 

(u2/f)'(x ) 

(u2/f)" (X ) 

(u2/f )'" {X ) 
IV 



( 



1 


A 


A 2 



\ 



(27) 



(u 2 /fY v (x ) 



\ (.n 2 /fp(x ) / 
where A n is the matrix from 



Art 





1 


A 




Proof. From theorem [2] we have that in (a, b) the functions u\ and u 2 defined by (jlip are 
solutions of the Sturm-Liouville equation P^|) satisfying the initial conditions (fT2"j) and (fTHj) , 
and the series in (jlip converge uniformly on [a, b]. Thus, the functions u\j f and u 2 /f are 
expended into uniformly convergent generalized Taylor series, and from (jllh we have that 

7fc(y)(a?o) = ~ 2 A2 



and 



,U2,, v l + (-l) fc+1 .ti 

7fc( j)(zo) = -j 1 A 2 . 
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Example 16 Our first example is the equation |1^[ ) with a constant coefficient q = —c 2 , c 6 
M., c / 0. Obviously, a nonvanishing solution of U0\) can be chosen in the form f{x) = e ca: . 
Considering xq = we calculate the first six rows of the matrix A n . We have 



A, 



On the other hand it is easy to see that the solution u\ from Ul\) in this case has the form 



1 1 



















1 
















-2c 


1 













4c 2 


-2c 


1 










-8c 3 


4c 2 


-4c 


1 




V o 


16c 4 


-8c 3 


12c 2 


-4c 1 


/ 



ui(x) 



C ~h K K ~. ^ 

-e KX + 



2k 



2k 



where k = \/ c 2 + A. Simple calculation involving the quotient u\j f = ^2k^ k c ^ x +%r e ^ ft+c ^ :r 
gives us the following values 



(«!//)' (0) = 0, (u 1 //) ,/ (0) = A, («!//)'" (0) 



-2cA, 



( Ul //) /y (0)=4c 2 A + A 2 , ( Ul /f) v (0) 



-8c 3 A - 4cA 2 



which is exactly the result of multiplication of the matrix A§ by the vector 

(1, 0, A, 0, A 2 , 0) T 

(see QT$). 

Example 17 Consider the equation 

u"{x) - a 2 x 2k u{x) = \u{x). 



(28) 



This equation with A = — 1 can be found, e.g., in Part 3, Chapter II, eq. 2.15] where it 
is stated that its solution is unknown, meanwhile the solution of i28\) when A = is known 



fWj Part 3, Chapter II, eq. 2.14] even for noninteger k. In particular, when k = — 2 2 " x and 
n is a negative integer number the general solution of the equation 



y"(x) - a 2 x 2k y(x) = 



(29) 



has the form 

y(x) = | x 2 "+! 



dx 



This solution can be used for solving 
general solution of A29\) takes the form 



c\ exp ( a(2n + l)x 2n + 1 J + C2 exp ( — a(2n + l)x 2n + 1 
Consider, for example, k = — 2 (n = 

a/x\ 



-I). The 



y(x) = ax(ae a l x 



c 2 e 
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As a particular solution f we can take, e.g., f(x) = axe a l x . Fix xq = 1. Calculation of the 
first five rows of the matrix A n gives us the following result 



( 1 




























V 




) 



2(a-l) 
8 - 16a + 4a 2 ^L e ~ 2a 1 

ft* 



1 



/ 



from where using theorem [75| the first four derivatives of the functions ui/f and u^j 'f at 
xq = 1 can be calculated. 

7 Conclusions 

A formula for calculating Taylor coefficients of the solution of the Sturm-Liouville equation 
is obtained. It is based on the knowledge of the Taylor coefficients of a particular solution 
of the equation corresponding to the zero value of a spectral parameter. The form of the 
(matrix) map transforming one set of Taylor coefficients into the other is independent of the 
particular form of the particular solution. 
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